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Abstract 

This paper addresses a detection problem where several spatially distributed sensors independently 
observe a time-inhomogeneous stochastic process. The task is to decide between two hypotheses regard- 
ing the statistics of the observed process at the end of a fixed time interval. In the proposed method, 
each of the sensors transmits once to a fusion center a locally processed sunmiary of its information in 
the form of a likelihood ratio. The fusion center then combines these messages to arrive at an optimal 
decision in the Neyman-Pearson framework. The approach is motivated by applications arising in the 
detection of mobile radioactive sources, and offers a pathway toward the development of novel fixed- 
interval detection algorithms that combine decentraUzed processing with optimal centralized decision 
making. 

Index Terms 

Decision making, sensor networks, inhomogeneous Poisson processes, nuclear detection. 

I. INTRODUCTION 

Decision making is crucial in translating information to action. Human decision makers can 
be significantly assisted in determining a time-critical plan of action if provided with concise. 
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dependable information. In a variety of applications — particularly those involving spatially and 
temporally nonuniform processes — the enhanced observational capabilities afforded by sensor 
networks render such solutions attractive for collecting the information required to arrive at a de- 
cision. Yet, the amount of data accumulated by a network of sensors is frequently overwhelming, 
underlining the need to filter and synthesize the collected data to ease decision making. Our goal 
in this work is to provide a framework that leverages local sensor-level information processing, 
to enable accurate fixed-interval decision making in the context of a spatially distributed sensor 
network which makes binary decisions based on observations of a time-inhomogeneous stochastic 
process. 

Sensor networks — both mobile and static — have been employed in a wide range of applica- 
tions, including environmental monitoring [[T|-[|3|, intruder detection [j4|, area coverage [[5|-[|9|, 



source localization p()|-p3|, and mapping of spatially distributed physical quantities p4| , p3|. 
Currently, networks of distributed sensors are used to trigger timely responses to a number 
of natural disasters, such as hurricanes p6| , earthquakes p7| and tsunamis [|T8|, p9|. In 
the classical approach to network-based decision making, sensors relay the entirety of their 
observations to a central processing unit, which analyzes the data and issues a global decision. 
While this centralized approach has the advantage of using all the information available, it 
does impose significant communication overhead. Alternatively, a decentralized decision-making 
scheme pO|-p3| can be used; in this setting, the sensors process measurement information and 
transmit a compressed version of it — typically in the form of a message with values in a finite 
alphabet — to a fusion center, which then provides a decision. Our motivation in this article 



stems from a class of problems associated with nuclear detection [24|-[27|; more specifically, 
the detection of illicit radioactive substances in transit. Remarkably, small or shielded quantities 
of nuclear material are very difficult to detect at a distance, due to the fact that their sensory 
signature is disguised in naturally occurring background radiation. Yet, the ability to provide 
fast and accurate decisions in such situations is of paramount importance to public safety and 



nuclear nonproliferation; see for instance [26|, [27|, which make explicit reference to the need 



for networks of detectors deployed along transportation routes. 



The physical quantities of interest in many applications (including nuclear detection [25|-[27|) 
can be captured by random processes characterized by discrete events that are highly localized in 
time. Phenomena of this sort can be mathematically modeled and analyzed within the framework 
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of point processes [28 l-BTJ. A realization of a point process is a random sequence of points, each 



representing the time and/or spatial location of an event. A point process can be characterized 
in terms of its intensity which corresponds to the rate at which events occur. Beyond nuclear 
measurement, typical examples of such processes include customers to and from a service facility 



in queueing theory [28 1, [|32|, electron emission from a photodetector in optical communications 



systems [33], generation of electrical pulses in neurons [34|, and others. Of special interest in 
nuclear detection are Poisson processes which provide the natural models describing the emission 
and measurement of radiation [|24), [|25J, [[35|. 

In regards to point processes, the problem of decision making between two alternative hy- 



potheses ("all clear" versus "alarm") has been addressed in [28 1, [36 1, [37 1; the solution typically 



involves the computation of a likelihood ratio, whose comparison against a threshold provides the 
decision. Error probability bounds for such decision problems are studied in [[38|, while robust 



decision making (in the presence of modeling uncertainty) is explored in [39|. Decision problems 



with time-inhomogeneous point processes also arise in optical communications [ |33| , [ [40| . In 
sum, for the classical (single observer) case, decision theory for general point processes is well- 
understood. However, the realization of these results in a network setting requires care. Indeed, 



the likelihood ratios in [28 1, [36 1 involve intensities computed on the basis of all accumulated 
information^ necessitating a modicum of caution in a setting such as ours where much of the 
computation and most of the raw data are decentralized (See Remark [2]). 

The problem of detecting (moving and stationary) radioactive sources using networks of 
sensors has received a fair bit of attention in the literature. In situations where the parameters 
(location, trajectory, activity) of the source are unknown, Bayesian methods are frequently used 
[ [T3| , [ [24| , [25[, [ |4T| , embedding the issue of detection in a parameter estimation problem. While 
powerful, Bayesian methods for source parameter estimation exhibit computational complexity 
exponential in the number of parameters estimated, posing challenges for their implementation 



in real time for networks with more than ten nodes [25 1, [41 1. An important insight — and one 
that serves as the starting point for our analysis — is that in many cases of interest, the problem 
of source localization can be decoupled from the problem of source detection. Indeed, there 



are improved methods [42[-[44[ for tracking the carrier of a potential radioactive source using 



Strictly speaking, the intensity is a conditional rate at which events occur (conditioned on available information). 
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sensor modalities other than a Geiger counter. Armed with this observation, source detection 
reduces to the problem of deciding whether the counts observed by a spatially distributed 
network of radiation sensors correspond solely to background radiation, or whether they also 



include emission from a radioactive source with known parameters. In this setting, [25| explores 
the Signal-to-Noise Ratio (SNR) resulting from the combination of data from a network of 
radiation sensors, allowing for spatially varying background rates. The analysis is restricted, 
however, to uniform linear source motion and does not provide a decision test. The costs and 
benefits of using networked sensors for moving sources, together with a threshold test (based on 
the total number of recorded counts) are addressed in [|45|, assuming uniform background and 
constant geometry between source and sensorj^For the case of a stationary source and correlated 
sensor measurements, a distributed detection scheme is developed in BS] using the theory of 



copulas. The work in [41 1 studies detection (via Bayesian estimation) for a moving source, but 



the motion is required to be linear with constant velocity. Detection and parameter estimation 



for an unknown number of static radioactive point sources are treated in p3| |, p4| . Evidently, 
the networked detection problem for general source motion with spatially varying background 
intensity has yet to be studied. 

Motivated by the above, we pose the following problem: a spatially distributed sensor network 
observes — over a fixed-time interval — a time-inhomogeneous point process which is known a 
priori to be governed by one of two intensities. How should the local information be pro- 
cessed and communicated through the network in order to reach a reliable decision regarding 
which intensity governs the observed process? With respect to the spectrum of approaches 
from centralized to decentralized, we take in this paper an intermediate approach that combines 
the significantly lower communication cost of decentralized processing (not decision) with the 
enhanced accuracy of centralized decision making. For the case of a vector of Poisson processes 
whose intensities explicitly depend on timej^ we develop an optimal — in the Neyman-Pearson 
sense — decision-making scheme that combines decentralized processing (local processing at each 
individual sensor) with centralized decision making via a fusion center. In particular, assuming 

^Our analysis indicates that the optimal test involves comparing the likelihood ratio against a threshold, rather than the total 
number of counts. 

^Time-dependence of intensities encodes the relative motion between the source and the sensors. 
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that the relative motion between the suspected source and the sensors is deterministic and known, 
and that the sensor observations are conditionally independent, our method relies on the sensors 
communicating processed information in the form of locally-computed likelihood ratios to the 
fusion center. The fusion center then combines these messages to arrive at a decision, without the 
need for any additional information such as the location or the raw data of individual sensors. 
As applied to radiation detection, our framework allows us to consider arbitrary continuous 
source motion in any number of dimensions allowing for sensor mobility and spatially varying 
background rates. In relation to sensor networks, the time-inhomogeneity in our problem leads 
to non-identically distributed sensor observations; this marks a departure from the frequently 
used independent identically distributed (i.i.d.) assumption. 

The paper is organized as follows. In Section |ll| we state the problem and our technical 
assumptions. The main result (Theorem [T]), which indicates how a global likelihood ratio test 
can be formulated based on local computation of sensor-specific likelihood ratios, is presented 



in Section III The proof of this result and the supporting technical material are found in Section 



IV Based on this analysis, we offer conservative lower and upper bounds on the probabilities 



of detection and false alarm, respectively, in Section IVl Finally, a numerical example of a 



one-dimensional case of networked nuclear detection is developed in Section VT highlighting 
the benefits of using multiple sensors. The results provided in this paper can be viewed as a 
building block toward a general decision-making framework that leverages networks of mobile 
sensor platforms to enhance detection capability in problems that involve time-inhomogeneous 
point processes. 

II. Problem Statement and Assumptions 

Consider a collection of sensors observing a point process generated by some physical phe- 
nomenon. The goal is to decide between two hypotheses regarding the state of the environment. 
To this end, each sensor communicates a processed version of its observations to a fusion 
center, which combines all received messages to a binary decision (Fig. [T]). With an eye towards 
applications such as detection of mobile radioactive sources where the point process can be 
observed only for a limited time, we require here that the decision be made within a fixed time 
interval. The problem is formulated as a binary hypothesis test based on measurements from an 
array of k sensors connected in a parallel network architecture. Figure [T] shows a realization of 
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such a network for radiation detection. 



moving source 




/ 

/ 

/ 



fusion center 

i 

decision 

Fig. 1. The architecture of the detection scheme as applied to radiation detection. There is a network of sensors each counting 
the total number of rays that have arrived at them. Thick arrows represent rays emitted from a moving radioactive source and 
thin arrows mark rays from background radiation. Based on the number and timing of arrival of those counts, each sensor 
computes a likelihood ratio L{i) which is then transmitted to a fusion center. The fusion center combines this information to 
make a decision regarding the presence of the radioactivity in the target moving in front of the sensors. 

Let Ho and Hi denote the two hypotheses regarding the state of the environment; for example, 
the absence or presence of a radioactive source on the moving target. 

Assumption 1: The reading at sensor i, where 1 < i < A;, is a non-decreasing Z+-valued 
piecewise constant, right-continuous function, which increases in steps of size one. 

Assumption 2: Conditioned on hypothesis Hj, j E {0, 1}, the observations at distinct sensors 
are independent. 

Sensor observations are modeled as inhomogeneous Poisson processes, whose intensities are 
deterministic functions of time. In particular, for I < i < k, sensor i observes a Poisson 
process, whose time-dependent intensity under Hq is while under Hi is + z/j(-). 

The explicit dependence of the intensities on time arises from the known motion of the target, 
which is a potential "source" of the point process with intensity Ui. Thus sensor observations 
are independent but not identically distributed. This is in agreement with the physics of the 
motivating application, since a gamma ray emitted from the source cannot pass through more 
than one sensor simultaneously (see Fig. [T]), and the time varying nature of the distance between 
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source and sensor changes the arrival statistics on the sensor side p5| . 
The following assumptions will be imposed on (3i and z/j. 

Assumption 3: For I < i < k, (3i : [0,T] ^ [/^min, /^max] is a bounded, continuous function 
with < /3min < /3max < oo, /3min, /^max independent of i e {1, 2, . . . , k}. 

Assumption 4: For 1 < i < k, Ui : [0,r] — )■ [z^min, ^^max] is a bounded, continuous function 
with < z/min < z^max < OO, z/^in, '^max independent of i G {1, 2, . . . , A;}. 



The physical significance of these assumptions will become clear in Section VI 
The detection problem can now be summarized as follows: Suppose that T > is the 
decision time; that is, the time by which a decision must be made. Given a single realization 
of a /c-dimensional vector of Poisson processes over the time horizon [0,T] (the k components 
corresponding to the k sensors), decide whether the intensities are given by the collection (3i(-) 
or by the collection + ^'i(-), I < i < k. 

111. Main Result 

In this section, we state our main result. Theorem [T] For a collection of sensors with a fusion 
center, configured in a parallel network architecture. Theorem [T] gives a procedure for locally 
processing sensor information and transmitting compressed summaries (at a single time) to the 
fusion center to enable networked decision making that recovers the performance of a centralized 
scheme. Our starting point is a measurable space (i7,^)j^on which a fc-dimensional vector of 
counting processes Nt = (A^i(l), . . . ,Nt{k)), t G [0,T] is defined. In our problem, Nt{i) is the 
number of counts registered at sensor zG{l,2,...,/c}upto (and including) time t G [0, T]. The 
two hypotheses Hq and Hi regarding the state of the environment correspond to two distinct 
probability measures on (fi,^). Hypothesis Hq corresponds to a probability measure Pq, with 
respect to which the Nt{i), I < i < k, are independent Poisson processes over t G [0,T] 
with intensities /3i{t), respectively. Hypothesis Hi corresponds to a probability measure Pi, with 
respect to which the Nt(i), I < i < k, are independent Poisson processes over t G [0,T] with 
intensities (3i{t) + z/j(t), respectively. The decision problem is thus one of identifying the correct 
probability measure (Pq versus Pi) on (Vl, ^) based on a realization of the fc-dimensional process 
iV, = (iV,(l),...,iV,(A;)). 

''Here, Q. is the sample space and ^ is a a-field on $7. Equipping (fi, ^) with a probability measure P gives the probability 
space (ri,^,P). 
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We will keep track of the flow of information using the filtration '■ < t <T) generated 
by the process Nt; here, for t G [0,T], = a{Ns : < s < t) is the smallest cr-field on Q 
with respect to which all the (/c-dimensional) random variables Ng, < s < t are measurable. 
The interpretation is: for any event A E , an observer of the sample path s A^^^ < s < t, 
knows at time t whether or not the event A has occurred. The a-field thus represents the 
information generated by the totality of sensor observations up to t = T; to wit, the information 
on which the decision must be based. 

A test for deciding between hypotheses Hq and Hi on the basis of observations can be 
thought of as a set Ai E with the following significance: if the outcome u E Ai, decide Hi, 
li oj E Aq = Vl \ Ai, decide Hq. For a test Ai E , two types of errors might occur. A "false 
alarm" occurs when the outcome u E Ai (i.e. decide Hi) while Hq is the correct hypothesis. 
A "miss" occurs when uj E Vl \ Ai (i.e. decide Hq) while Hi is the correct hypothesis. Clearly, 
the probability of false alarm is given by Po(Ai), while the probability of a miss is given by 
Pi(r] \ v4i). Then, the probability of detection is given by Pi(v4i) = 1 - Pi(fi \ Ai). 

In the Neyman-Pearson framework, one is given an acceptable upper bound on the probability 
of false alarm a E (0, 1), and the problem is to find an optimal testj^ a set A\ E which 
maximizes the probability of detection over all tests whose probability of false alarm is less 
than or equal to a. The following result provides an optimal test that employs local information 
processing at the sensor level, to enable decisions at the fusion center that recover the optimal 
performance of a centralized Neyman-Pearson test. 

Theorem 1 (Main Result): Consider a network with k sensors and a fusion center connected 
in the parallel configuration of Fig. [Tj For 1 < i < /c, let Nt{i), t E [0, T] denote the observation 
at sensor i over the time interval [0,T] and let (r„(z) : n > 1) be the jump times of Nt{i). 
Assume that at decision time T, sensor i transmits to the fusion center the statistic 

M.)^e.p(-f.,M.)|'(l.|(||) 

computed on the basis of its observation t ^ Nt{i), t E [0,T]. Then, the test A\ = {Lt > 7} 



We restrict attention iiere to tests without randomization; see 



46 
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performed at the fusion center, with 



k 



i=l 



and 7 > satisfying Fq{Lt > 7) = aj^ is optimal for ^^-observations in the sense that for 
any Ai G with Po(Ai) < a, we have Fi{Al) > Pi(Ai). 

Before continuing with the proof of Theorem [1] we highlight why optimal decision making 
through decentralized processing of information at the sensor level is possible in the case 
considered here. 

Remark 1: The decision test {L^ > 7} is optimal (in the Neyman-Pearson sense) for 
observations, the latter comprising the totality of information in the k waveforms t Nt{i), 
t G [0,T]. Equivalently, if one were to consider a centralized framework (where information is 
continuously streamed from each of the sensors to the fusion center), {Lt > 7} would be the 
optimal test. Note, however, that each Lr(i) in the product Lt = Y[^=i^T{'i) can be computed 
locally at sensor i without any knowledge of the measurements at other sensors (see also Remark 
[2]). Indeed, computation of Lt^i) requires knowledge solely of the times (r„(i) : n > 1) at which 
counts have been recorded at sensor i, and the quantities /3i(t), i'i{t), which are deterministic. 
Consequently, each sensor simply needs to transmit its locally computed Lt^i) to the fusion 
center at (the single) time t = T (in lieu of the element t Nt{i), t E [0, T] of function space). 
The fusion center then forms the product Lt = Y[i=i -^t(^) which is compared against 7 to arrive 
at a decision. We thus retain the accuracy of centralized decision making while decentralizing 
most of the data processing, thereby accruing significant savings in communication costs. 

IV. Proof of Main Result 



The contents of this section are organized as follows. In Section IV-A we provide precise 



definitions for various quantities of interest. In Section IV-B, we state and prove some results 



needed for the proof of Theorem [T] In Section IV-C the proof of Theorem [T] is completed 



^While there may not exist such 7 for every a G (0, 1), one can always find a sequence a„ — >■ for which there exist 7,1 > 
with ¥o{Lt > 7n) = Qn- In other words, arbitrarily small upper bounds on probability of false alarm can be accommodated. 
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A. Definitions 

Let ^,P) be a probability space. The sample space Vt is the set of all possible outcomes 
w of a random experiment, ^ is the cr-field (or cr-algebra) of events, and P is a probability 
measure. A point process on [0, oo) can be described by a sequence (r„ : n > 0) of random 
variables defined on (r2,^,P) taking values in [0,oo] such that 

To = , (1) 

r„ < OO =^ r„ < r„+i . (2) 

Here, t;„ denotes the (random) time of the n-th occurrence of an event (such as a radiation 
counter registering a count). Associated to the sequence (r„ : n > 0) is the stochastic process 
{Nt:t> 0) defined by 

Nt = Yl Mr.<t) , (3) 
re>l 

where 1a denotes the indicator function of A, i.e. 



1 if CO e A, 
else. 

Thus, Nt counts the number of occurrences of the phenomenon prior to or at time t. Nt is a 
counting process, i.e. Nt is a Z+-valued process with A'^o = such that the sample paths t Nt 
are non-decreasing, piecewise constant, right-continuous functions of t which increase in steps 



of size 1. We will also refer to Nt as a point process [28 1. 

A filtration {^t : t > 0) is an increasing family of sub-a-fields of i.e. C ^ for all 
t, and s < t implies C ^t- The cr-field ^t represents the information available at time t. 
A stochastic process {Nt : t > 0) taking values in M'"' is said to be adapted to the filtration ^t 
if for all t > Nt is -measurable; i.e. for any Borel measurable subset B of M^, the event 
{Nt e B} e ^ox t > 0, let = a(iV, : < s < t) be the smallest cr-field on Vt with 
respect to which all the random variables A^^, < s < t are measurable. Then, {^^^ : ^ > 0) 
is the filtration generated by the process Nt and corresponds to the information available to an 
observer of the process Nt. Clearly, if Nt is =^t-adapted, then C ^t for all t. Next, let us 
make precise what we mean by an inhomogeneous Poisson process. 
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Definition 1 (Inhomogeneous Poisson process): Suppose \{t) is a nonnegative, measurable 
function such that j^\{s)ds < oo for all t > 0. A point process Nt on a probability space 
(fi,^,P) adapted to the filtration : t > 0) is said to be a (P, #j)-Poisson process with 
intensity \{t) if for < s < t, 

1) Nt — Ngh independent of and 

2) Nt — A^s is a Poisson random variable with parameter J^* \{u)du, i.e. for all n G Z+, 



/*A(r)dr 

¥{Nt -Ns=n) = e-^s Hr)dr\^ ^ /_ _ 



B. Useful Results 

The primary result in this section is Proposition |2} which provides the probabilistic setup for 
the statement and proof (given in the next section) of Theorem [1} Our development proceeds 
through the following steps. We start with Proposition [l] — a result of Bremaud [28| which plays 



a fundamental role in our analysis. To apply Proposition [T| we first relate our problem — posed on 
the time interval [0, T] — to a corresponding problem over the time interval [0, 1], as in Proposition 
[T| This is accomplished by a time-rescaling argument (see proof of Proposition |2]). Lemma [2] 
aids us in this regard by describing how the intensity of a Poisson process transforms under a 
change of time. Next, we verify in Lemma [1] that condition ([8]) in the statement of Proposition 
[T] holds. It is important to note that Proposition [T] and Lemmas [2] and |3] are formulated in 
terms of a general probability space, not necessarily identical to the one in Proposition |2] (which 
supports our processes of interest). Lemma [1] on the other hand, pertains to the specific setup 
of Proposition |2j Finally, Lemma |3] recalls the Neyman-Pearson Lemma, which is used in the 
next section in the proof of Theorem [T] 

In allowing stochastic intensities, as in Proposition [T] below, one has to have the technical 
requirement of predictability [ |28| Section 1.3]. In our problem, however, the intensities are 
deterministic and automatically predictable. In terms of notation, we follow the convention of 
using \t for a stochastic intensity versus \(t) for a deterministic one. 

Proposition 1 (Theorem VI.2.T3, 128]): Let (Nt{l) , . . . , Nt{k)) be a fc-variate point process 
adapted to the filtration on the given probability space (fi,^,Po), and let Xtii), I < i < k, 
be the predictable (Pq, ^t) -intensities of Nt{i), I < i < k, respectively. Let I < i < k. 
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be nonnegative, -predictable processes such that for all t > 0, 1 < i < A;, 

fis{i) Xs{'i)ds < oo Fq — a.s. (5) 
Let {Tn{i) : n > 1) denote the jump times of Nt{i), I < i < k, and define the process Lt by 

k 

L, = , (6) 



i=l 



where 




^tW =111 /^r„(i)W I ^^p{^ - ^s{i))\s{i)ds^ , (7) 

with the convention that nn=i(- • • ) = 1- Suppose moreover that 

Eo[Li] = 1 . (8) 
Define the probability measure Pi on (Q,^) by 

Then, for each I < i < k, Nt{i) has the (Pi, -intensity \t{i) = iJit{i)\{i) over [0,1]. 

To underscore the role that deterministic intensities play in allowing decentralized processing 
in our problem, we briefly discuss some of the subtleties that arise in networked detection of 
point processes with stochastic intensities. 

Remark 2: In Proposition [Tj the assumption of -predictability of \t{i), i^tii) (which are 
in general stochastic) implies that the latter potentially depend on all the information in ^t, 
which includes the information generated by all the sample paths s ^-^ Ns{i), s E [0,t], 1 < 
i < k^ Consequently, for a sensor network which observes a point process with stochastic 
intensities, computation of the Lt(i)'s cannot be decentralized as described in Remark [T] without 



incorporating a filtering component |28 Section VI.4] — finding the best estimates Lt{i) of Lt{i) 
based on the locally available information ^^^^^ = a{Ns{i) : < s < t). In such a decentralized 
processing scheme, the likelihood ratio to be compared to a threshold (at the fusion center) would 
be the process Lt = Y['i=i^t{i) evaluated at the decision time. This would entail, in general, 
some loss of performance in comparison to a fully centralized scheme. These extra considerations 

may contain information generated by other random quantities too. 
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do not arise in our detection problem since our analysis focuses on a problem with deterministic 
intensities. 

Next, we state Lemmas [T] and [2} which together verify ([8]) for our problem. The proofs of 
these Lemmas are given in the Appendix. Lemma [T] shows that the process (Lt : t E [0,T]) 
defined by ([TOl)-([TT]) is a martingale with mean one. Lemma |2] describes how the intensity of a 
Poisson process transforms under a time rescaling. Taken together, these lemmas enable us in 
Proposition |2] to obtain the probability measure Pi on ^) (corresponding to hypothesis Hi). 
In the sequel, we denote the expectations corresponding to probability measures Pq and Pi by 
Eo and Ei, respectively. 

Lemma 1: Under Assumptions Isl |4l the process (Lt : t E [0,T]) defined by 



k 

Lt^l[LS) (10) 

i=l 

where 



n=l 

is a nonnegative (Pq, =^/^)-martingale. By (Pq, ^j^) -martingale, we mean that 

1) Lj is adapted to the filtration 

2) Eo[\Lt\] < oo for all t E [0,T], 

3) For < s < t < T, Eo[Lt|^f ] = L„ Po-a.s. 

Thus, Lt has constant mean, i.e. Eo[Lf] = Eo[Lo] = 1 for all t E [0,T]. 

Proof: See Appendix. ■ 

Suppose (i7,=^,P) is a probability space equipped with the filtration : t >0). Let A(t) 
be a nonnegative, measurable function defined on [0, oo) with \(s)ds < oo for all t > 0. Let 
{Xt : t > 0) be a (P, ^() -Poisson process with intensity X{t) (see Definition [TJ. 

Lemma 2: Fix T > 0. Let m = t/T. Let = Xt = Xt-u, ^u = '^t = ^t-u for m > 0. Define 
\{u) on [0, oo) by 

\{u)=T\{J-u) (12) 

for M > 0. Then, (F„ : m > 0) is a (P, ^„)-Poisson process with intensity \{u). 

Proof: See Appendix. ■ 
We now state Proposition |2} which constructs a probability measure Pi on (i^,=^) corre- 
sponding to hypothesis H\. This completes the construction of the probabilistic model of our 
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networked decision problem. As will be seen in the next section, the explicit construction of Pi 
via the process Lt and the probability measure Pq, as described in the proof of Proposition |2} 
facilitates the application of Lemma [3] in proving Theorem [T| 

Proposition 2: Suppose (fi, =^,Po) is a probability space, on which Nt = (A^t(l), • • • , Nt{k)), 
t E [0,T], is a vector of independent ^/^-Poisson processes whose components Nt{i) admit 
intensities (3i{t), 1 < i < k, with (3i{t) satisfying Assumption |3] Then, there exists a probability 
measure Pi on (i^, with Pi <^ Pq, with respect to which Nt{i), I < i < k, are indepen- 
dent ^/^-Poisson processes over t E [0,T] with intensities l3i{t) + z/j(it), with z/j(t) satisfying 
Assumption |4} 

Proof: Recall the process {Lt : t E [0,T]) defined by ([10]), ( [TT] ). Lemma [T] assures us that 



Lt is a nonnegative random variable with EofL^] = 1. Hence, Pi defined through 

Pi(A) = [ LT{uj)Fo{dcj) for Ae,^, 

J A 

is indeed a probability measure on {Vl, ^) which is absolutely continuous with respect to Pq. We 
would now like to show that with respect to Pi, Nt{i) iox 1 < i < k, are independent Poisson 
processes over t E [0,T] with intensities (3i{t) + h'i{t). To enable the application of Proposition 
[T] we use a time rescaling argument. 

Let u = t/T be a rescaled time variable taking values in [0, 1]. For 1 < i < k, u E [Q, 1], 
let iV„(«) = Nt{T) = NT.u{^). Let iV„ = (iV„(l), . . . , iV„(A;)) and let ^„ ^ cr{N, : < v < u) 
for u E [0, 1]. By Lemma [2} each Nu{i), 1 < i < k is a (Pq, ^M)-Poisson process with intensity 
I3i{u) = T(3i{T ■ u). The independence of iV„(i) — iV„(i) and iVp(j) — Ng(j), i ^ j follows from 
the independence of Nx-uii) — Nx-vii) and Nt-pU) — Nr-qij). Denote by (f„(i) : n > 1) the 
sequence of jump times of Nu{i). Note that f„(z) = Tn{i)/T. Letting Vii^u) = Tui(T ■ u) for 
I <i <k, define a process [L^ : u E [0, 1]) by 

k 

Lu = \[L^{i) (13) 

i=l 

with 

^ exp (- r i.,{s)ds) (1 + • (14) 

\ Jo J „=i V Pi[Tn{l))J 

It is now easily checked that Lt = Lt/x for all t E [0,T], which implies in particular that 

^=i. ■ (15) 

October 5, 2012 DRAFT 



C. D. PAHLAJANI, I. POULAKAKIS, H. G. TANNER 15 

We now apply Proposition [T] [[28| Theorem VL2.T3], using the rescaled time variable u E [0, 1], 
filtration with $i{u) and l + z/j(u)//3j(n) in place of \t{i) and jitii) respectively, to infer that 
Nu{i) has intensity f3i{u) + with respect to Pi. Using Lemma [2] in the "reverse" direction 
(i.e. interchanging t and u, replacing T by 1/T), it now follows that the A^t(i)'s, I < i < k, 
have (Pi, ^/^) -intensities (3i{t) + z/j(t), respectively. To complete the proof, it remains to show 



that the Nt{i), 1 < i < k, are. independent and Poisson under Pi. Using Assertion of [28 
Theorem II.3.T8] (with Xt = 1), we get that for 1 < i < A;, 



Nt{t)- [ [p,{s) + u,{s)]ds (16) 
Jo 



is a (Pi, ^/^)-martingale. By the Multichannel Watanabe Theorem [28 Theorem II.2.T6], it 
now follows that with respect to Pi, the A^t(i)'s are independent ^/^-Poisson processes over 
t E [0,T] with intensities /3j(t) + Ui^t), respectively. ■ 

Before concluding this section, we state the Neyman-Pearson Lemma which will be used in 
the next section to prove Theorem [Tj The Neyman-Pearson Lemma describes an optimal rule for 
deciding between probability measures Pq and Pi on a measurable space (Vl, ^) on the basis 
of observations in the sub-a-field ^ C . Thus, for any event A G ^, it is known whether or 
not A has occurred. Recall that Pi is absolutely continuous with respect to Pq, both restricted 
to (^],^), denoted Pi < Pq if, whenever Ae^ with Po(A) = 0, we have Pi(A) = 0. 

Lemma 3 (Neyman-Pearson Lemma, Theorem VLLTl, /[2^j." For a E (0, 1), suppose 7 is a 
real number such that 

Po(L > 7) = « (IV) 

where L is the Radon-Nikodym derivative of Pi with respect to Pq, both probabilities on [Vt, 
Then the decision strategy A\ = {L > 7} is optimal for ^-observations in the sense that for 
any Ai G ^ with Po(Ai) < a, we have Pi(A^) > Pi(Ai). 

See [ |28| for a proof. We will use this Lemma in the next section to prove Theorem [T] with 
playing the role of ^. 

C. Proof of Theorem [7] 

Let Po, Pi be the restrictions of Pq, Pi respectively to . Since Pi < Pq on 3^, the 
restrictions of Pi and Pq to the smaller a-field inherit the absolute continuity, i.e. Pi ^ Pq. 
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Hence, the Radon-Nikodym derivative exists; i.e. there exists a nonnegative, =^^-measurable 
random variable ^, denoted such that for any A G 



F^{A) = / a^)Foidu) . (18) 

J A 

Moreover, this Radon-Nikodym derivative is unique in the sense that if ^ is any nonnegative. 



-measurable random variable satisfying (18) with ^ replacing ^, then ^ = ^, Po-a.s. Since 
Lt is a nonnegative, ^^-measurable random variable (by Lemma [T]) which satisfies (18), it 
follows that 

Po-a.s. A direct application of Lemma [3] completes the proof. 

V. Performance Analysis 

Here we provide a lower bound on the probability of detection and an upper bound on the 
probability of false alarm when the proposed detection scheme is used. It turns out that bounds 
on both these probabilities involve the tails of (different) Poisson distributions. To compactly 
describe our results, we follow the notation of [48 1. 

Definition 2 (Poisson Tails): For A > 0, j G Z+, let p(A, j) denote the Poisson distribution 

P(A,J) = e"^^ • 
The left and right tail probabilities are defined by 

n oo 

P(A,n) = Y,p{X,j) , P(A,n) = ^p(A,j) , (20) 

j=0 j=n 

respectively. Note that P(\,n — 1) + P{\,n) = 1. 
The following quantities will also be of interest: 

k „T k „T 

B = y2 '^*(^)^^ ' J = 11 ''^(^)^^ • (21) 

i=l ^0 Jo 



It now follows from ([TO])-([TT) that 

k Nrii) 



i=l n=l 

In the sequel we will also use the integer ceiling function [■] which assigns to a real number x 
the smallest integer greater than or equal to x. 
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For 7 > 0, consider the test A* = {Lt > 7}- A lower bound on the probability of detection 
IPi(-^T > 7), and an upper bound on the probability of false alarm Fo{Lt > 7) can now be 
obtained as follows. Recalling Assumptions |3] and |4} define 



C = l + 
D = l + 



< min min 1 + 

l<i<fc 0<t<T 



m 

> max max 1 + — — - 

/3min ^<i<k 0<t<T \ (3i{t) 



(22a) 
(22b) 



Note that if one can find such that i~ < Lt < i^, then for j G {0, 1}, 

Pi(^" > 7) < Pj(^T > 7) < Pi(^^ > 7) • 

Letting 



it can be verified that 



i=l 

i- <Lt< 



Nrii) 



Nrii) 



i=l 



. Next, note that 

k 

>7 ^ E^^(^) 



> 



i=l 
k 



>7 



l) > 



i=l 



log 7 + J 

log 7 + J 
logD 



Since A^t(^) for 1 < i < k are independent Poisson random variables with parameters /3i{s)ds 
with respect to Pq, it follows that Yl^=i ^T^i) is a Poisson random variable with parameter B = 
Yl'i=i lo l^i{^)ds with respect to Pq. Under the probability measure Pi, A^r(^) for 1 < i < A; are 
independent Poisson random variables with parameters /Q^[/3j(s) + Vi{s)\ds. It follows that under 
Pi, XlLi ^T{i) is a Poisson random variable with parameter J + B = Xli=i Jq [i^iis) + /3i{s)]ds. 
Hence, 

"log7 + J 



Pi(^T >i)>pyJ + B 
MLt >i)<p (b, 



logC 
log 7 + J 
logD 



VI. Application to Nuclear Detection 



(23a) 
(23b) 



In this section, the framework developed above is applied to the problem of detecting ra- 
dioactive materials in transit, using a network of spatially distributed sensors. The setting here 
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is simple but representative of a frequently encountered class of scenarios. Our method is not 
restricted, however, to this setting. Indeed, the results in the previous sections apply whenever 
the intensity of the suspected source and the motion of the source relative to the sensors are 
deterministic and known. 

Radiation sensors always record background radiation (due to cosmic radiation and due to 
naturally occurring radioactive isotopes in the environment). In the absence of illicit nuclear 
material (hypothesis Hq is true), the sensors simply measure background. If radioactive material 
is present (hypothesis Hi is true), the sensors record the sum of the photons coming from 
background and the photons coming from the material. These two sources of radiation act 
independently, and one can treat each sensor as observing a single Poisson process whose 
intensity is the sum of intensities due to background and material (the source). The problem 
we face is to determine, in a fixed amount of time, whether a target passing in front of the 
sensors is a source of radiation. 

The specific assumptions for this problem are as follows: The workspace is the horizontal 
plane, M^. We have k sensors uniformly spaced along the positive x-axis at locations x = 0, 
X = i, X = 2£, . . . , X = (k — l)i in a configuration as that shown in Fig. [T| To span a length 
of approximately 100 m, we choose i = 11 m, k = 10, and for simplicity, we assume that 
the sensors are identical. Let t E [0, T] denote time, where t = corresponds to the instant 
the count recording is initiated, and t = T is the final time at which a decision regarding the 
existence of a source is to be made. Let /3j(t) be the intensity of background radiation at the 
location of sensor i, 1 < i < k, which does not have to be uniform and in general can be time- 
dependent. For simplicity, we assume in this example that background intensity is time-invariant, 
so /3j(t) = /3j G N, where (3i is assumed to be varying between locations, from a minimum of 
/3min = 2 counts per second (cps) to a maximum of /3max = 8 cps, with the maximum appearing 



at the first and last sensor and the minimum occurring at the sensor in the middle (Fig. 2(a)). 
We assume that a target is passing at a distance h = 0.362 m (the equivalent of 14^ inches) 
from the x-axis, namely with a constant coordinate y = h, appearing first at some initial location 
(xo, h) = (—4, 0.362) G M^, and moving with constant speed v = 17 m/s (roughly 38 mph) in 
the direction of the positive x-axis. 

To illustrate the derivation process, let us for the sake of argument assume that the acceptable 



probability of false alarm in this scenario is a = 10 (see ( |23b[ )). With rj(t), 1 < i < k. 
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denoting the distance between sensor i and the potential source, the intensity z/j(t) at sensor i 
due to the source is modeled in \25j b>j^ 



_Xa_ 



(24) 



where a > is the activity of the potential source (in cps) and % > (in m^) is the sensors' 
cross-section coefficient]^ We assume a numerical value for equal to what has been used in 



|25|, but shielded in 3 cm of lead, dropping the source's perceived intensity by one order of 
magnitude to = 506.8 cps-m^. We also assume that no sensor is ever closer than distance h 
to the target, ensuring that Vi(t) is always bounded. 

Since the location of the potential source at time t E [0, T] is {xo + vt, h), the distance rj(t) be- 
tween the potential source and sensor i, 1 < i < k is given by rj(t) = a/ {xq + vt — {i — + h'^. 
Recalling (|24]) with 

^ (k- l)i - 2xo 

T = ^ ^ = 6.3 s 



for the decision time, we get 
for t e [0,T]. Since 



{xo + vt-{i- i)ey + 



T 



Vi{s)ds 



from (21 1 we obtain 



J 



hv 



hv 



tan 



xo + vT-{i- l)i 



i=l 



tan 



h 



-1 I xo+vT-{i- 
h 



tan 



1 f X q - (i - l)i 
h 



and for the ten-sensor array we have J = 2559.74 counts. 
It is not hard to see that 



Tmin = min min rAt) = h 

l<i<k 0<t<T 



fmax = max max 

l<i<k 0<t<T 



^In fact, for a planar detection scenario the sold angle scales proportionally to l/r^. 
'For the case of heterogenous sensors, each will have its own Xi- 
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Fig. 2. Distribution of background intensity \2(a)\ and integrated perceived source intensity ( |2(b)[ ( over the considered ten-sensor 
array. 



Thus, for the case of the ten-sensor array, ([22]) evaluates to 



c 

D 



[{xo + vTy + h^]f3, 



5.97 X 10" 



1 + 935.24 . 



With reference to (21 1 and Fig. 2(a) 



we have B = counts, and with this we can attempt to 



numerically compute a threshold 7 for the likelihood ratio test using (23b). It can be verified that 



the Poisson tail on the right hand side of ( |23b[ ) falls below 10 when the second argument of 

^ >338 



P(-, •) increases to 338 (see Fig. 3(a) ). We thus compute the value of 7 for which 
counts, and obtain that with 7 = 0.1718, the bound on the probability of false alarm Pq falls at 
8.5 X 10^^, which is below the acceptable error rate. The decision rule therefore is based on the 
test: 

Li > 0.1718 (25) 
which if true, suggests that the target is indeed a radioactive source. It should be mentioned that 



there is conservatism built in the bounds (22), which renders the probability of detection using 



( |25] ) rather impractically small for the given false alarm rate. In addition, it is acknowledged 
that the illustrated method for obtaining a threshold makes the solution for 7 very sensitive 



to changes in the underlying parameters h, v, x*^ and /Smin- Improving the bounds in ( [22] ) is 
part of ongoing work. Nevertheless, the analysis still gives insight into the effect of different 
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(a) Poisson tail 



(b) False alarm rate 



Fig. 3. Bound on the probability of false alarm (PFA). The Poisson tail that upper bounds the probability of false alarm 
decreases monotonically with its second argument (Fig. |3(a)^ , and a value for the latter can be identified for which the former 
falls below a desired value. Without resetting the threshold constant, we see that the upper bound on the PFA decreases with 
the addition of new sensors. For fc = 8 sensors and 7 — 0.1718 a PFA is at most 20%, while for two additional sensors one 
can guarantee a one-in-a-million chance for a false alarm. 



parameters on the probability of detection. To illustrate that point, let us consider the possibility 
of using more sensors with the same spacing £ as before. Without changing the decision rule 
( |25| ) (keeping the same threshold), the analysis shows (Fig. |3(b)[ ) how the upper bound on the 



probability of false alarm (PFA) estimated in (23b) not only falls monotonically with the addition 



of new sensors, but that there is a clear transition between the state where the sensor network 
decision is unreliable, and that where an alarm should be taken into account seriously. Such 
information can be useful for determining the maximum number of detectors that can get out 
of commission before significantly compromising the effectiveness of the system. 

VII. Conclusions 

A network of sensors can be deployed to optimally decide between two hypotheses regarding 
the statistics of a time-inhomogeneous point process in a way that preserves the accuracy of 
centralized decision making without incurring the increased communication cost. The sensors 
collect their measurements over a fixed-time interval, at the end of which a processed summary 
is communicated to a fusion center. In particular, each sensor transmits a locally computed 
likelihood ratio to the fusion center, which then compares the product of the sensor-specific 
hkelihood ratios against a threshold to arrive at a decision. The analysis is based on the Neyman- 
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Pearson formulation. A set of conservative performance bounds on the error probabilities is 
provided and the framework is applied to the problem of detecting a moving radioactive source 
using an array of sensors. The work here supports the development of a general decision- making 
framework that leverages networks of mobile sensor platforms to enhance detection capability 
in problems that involve time-inhomogeneous point processes. 

VIII. Appendix 

Proof of Lemma [7]- The process Lt given by pO] ) and ( [TT] ) is a particularization to our 
problem of the general process Lt in Proposition [T] The latter admits the representation [28 
Equation VI. 2. 4] 

[i] - l)dM,{i) (26) 



where Mt{i) = Nt{i) — \s{i)ds, and for t > 0, Lt- = limgy^tLs is the left limit of Lt. The 



application of (26) to our problem, with Lt given by ( [T0| ), ( [TT] ), yields 



i=l 







(27) 



where Mt{i) = Nt{i) — L f3i{s)ds. The non-negativity of Lt is evident from ( [10] ), ( [TT] ). To 



complete the proof, it thus suffices to show that each of the integrals on the right in ( [27] ) is a 

ft 



martingale. By [|28| Theorem II.3.T8], for 1 < i < k, Ls-{Pi{s) / Pi{s))dMs{i) is a (Po,=^ 
martingale whenever 



En 



Ls- Ui{s)ds 



< oo 



(28) 



for t > 0{^By Assumptions [sj [ij we get that 

k 

Lt<l[K 

i=l 



Ntii) 



(29) 



28 



'"Actually, Theorem 1I.3.T8 in 
■^/'-adaptedness of Lt-, by 28 Theorem I.3.T5] 



also requires that Lt- be -predictable. This follows from the left-continuity and 
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where K = 1 + z/max//3min- Since K > 1 and the Nt{i), 1 < i < k are independent with each 
Nt{i) non-decreasing in t, we get 



En 



Ls- Ui{s)ds 



i=l 
k 

< I^max t J]^ exp 

i=l 



{K-l) / Pi{s)ds 



< oo 



(30) 



(31) 



where the last line follows from the fact that under Pq, for t > 0, each Nt{i) is a Poisson random 
variable with parameter (3i{s)ds. ■ 
Proof of Lemma 2- Note that \{u) is nonnegative and measurable with \{v)dv < oo for 
all M > 0. Next, since Xt-u is =^T M-measurable for all m > 0, it follows that Yu is ^^-adapted. 
To complete the proof, we need to show that iox < v < u, — Y^ is independent of and 
is a Poisson random variable with parameter \{r)dT. Since Xt-u — Xt-v is independent of 
^T-v, it follows that y„ — y„ is independent of W^. Finally, for n G Z+, we have 



P(F„-n = n) = F{Xt.u - Xt.. = n] 



exp 



exp 



\{s)ds 



T-v 



ff'^\{s)d^ 



\{T)dT 



i:kr)d 



where the last equality follows by making the change of variables r = s/T. 
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